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Abstract 

We extend our previous studies of the neutrino-induced reactions on 12 C 
and muon capture to include partial occupation of nuclear subshells in the 
framework of the continuum random phase approximation. We find, in con- 
trast to the work by Auerbach et al., that a partial occupation of the pi/ 2 
subshell reduces the inclusive cross sections only slightly. The extended 
model describes the muon capture rate and the 12 C(u e ,e~) 12 N cross sec- 
tion very well. The recently updated flux and the improved model bring 
the calculated 12 C(f^, /x~) 12 N cross section (~ 17.5 • 10 cm 2 ) and the data 
((12.4±0.3(stat.)±1.8(syst.)-10 -40 cm 2 ) closer together, but does not remove 
the discrepancy fully. 

25.30.Pt 



1 



Typeset using REVIgX 



I. INTRODUCTION 



In recent years the study of neutrino-nucleus reactions has attracted continuously growing 
interest. Besides astrophysical motivations this interest has particularly been inspired by 
the use of neutrino-nucleus reactions as a laboratory to study fundamental properties of 
neutrinos. The best studied nucleus by z^-induced reactions is 12 C, which is an ingredient of 
liquid scintillator and plays an important role in the oscillation search experiments performed 
by the KARMEN § and LSND § collaborations. 

The exclusive and inclusive charged- current reactions l2 C(v e , e~) 12 N have been measured 
by the KARMEN and LSND collaborations for neutrinos generated by muon decay at rest 
The KARMEN collaboration has supplemented this by a measurement of the neutral 
current 12 C(z/, z/) 12 C*(15.11 MeV) cross section [Q. Using the LAMPF pion-decay-in-flight 
neutrino source, the LSND collaboration has also measured the 12 C(z^, //~) 12 N cross section 
||. Together with the precisely known inclusive and exclusive muon-capture rates on 12 C 
0, these measurements represent a challenge to the theorists. 

Various data on the exclusive reactions (also including electro-magnetic transitions and 
beta decays), leading to the triad of T = 1, J n = 1 + analog states in the A = 12 nuclei 
(these are the 12 B and 12 N ground states and the excited state in 12 C at 15.11 MeV), 
appear to be consistent and are theoretically well described, validating the total neutrino flux 
determination ||. However, the situation is somewhat different for the inclusive reactions; 
the term is adopted here for the reactions leading to final states in 12 B and 12 N other than the 
ground state (in the following denoted by 12 B* and 12 N*, respectively) and involving states 
in the continuum. The data on inclusive weak processes ( 12 C(^ e , e~) 12 N*, 12 C(f M , /i~) 12 N* 
and muon capture) are usually not simultaneously reproduced by calculations @,|Hj, with 
a notable exception (Tl|. For example, a reasonable model to calculate these cross sections 
is the continuum random phase approximation (RPA). In fact, one finds that this model 
describes the inclusive muon capture rate and the 12 C(u e , e~) 12 N* cross section very well, but 
it overestimates the 12 C(z/ A1 , /i~) 12 N* data noticeably. Of course, these three weak processes 
probe the nuclear response at different momentum (q) and energy (u) transfers, and thus 
the good description of the capture rate (with average values To ~ 25 MeV, q ~ 80 MeV/c) 
and the f e -induced cross section (uJ « 23 MeV, q ~ 50 MeV/c) might simply indicate that 
the model does well for low values of q and u, but it is less reliable at larger momentum 
and energy transfers which are important for the l2 C(h> fl , /i~) 12 N* experiment (cZJ ~ 37 MeV, 
q « 200 MeV). Such hypothesis has, however, been invalidated by the continuum RPA 
study of the inelastic electron scattering on 12 C where the model has been demonstrated 
to describe well the data for momentum and energy transfer regions relevant for the LSND 
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Could the origin of the discrepancy be nuclear correlations missing in the continuum RPA 



model? This is in fact what the authors of Ref. |TT| suggest. First, they point out that, 
due to pairing effects, there are about 1.6 nucleons occupying the pi/2 subshell in the 12 C 
ground state, rather than the closed p 3 / 2 configuration (and empty pi/ 2 subshell) assumed 
in the previous continuum RPA calculations. Second, they argue on the basis of a standard 
RPA calculation that this partial occupation in fact allows one to simultaneously describe 
all three inclusive weak processes. 

Obviously the method of choice to test this conjecture is the interacting shell model. 
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But a reliable calculation of the 12 C(i/ M , /i~) 12 N* cross section requires inclusion of at least 
the complete 3hu model space, which is a formidable task even on modern computers. 
Short of such a large-scale shell model calculation, we have performed a set of continuum 
(standard) RPA calculations for the various inclusive (and exclusive) cross sections extending 
our previous model to allow for partial occupancies of the subshells. The second motivation 
for repeating our previous calculations within an improved model is that recently the decay- 
in-flight neutrino spectrum for the LSND experiment has been revised. With the updated 
flux, the experimental cross section slightly increased to (12.4 ± 0.3(stat.) ± 1.8(syst.) • 10 -40 
cm 2 ]T3| . As the LSND experiment is such an important benchmark we feel obliged to also 



consistently update our calculated cross section allowing for a meaningful comparison. 

The continuum random phase approximation for nuclei with closed-shell configurations 
is well documented in the literature [[L4[]. A generalization of this formalism to partial 



occupancies is readily achieved by multiplying the relevant matrix elements in Eq. (17) of 
Ref. |TJJ with the partial occupation numbers, rih, for the hole states, thus 



(ph^iyip'h'- 1 ) - ^^(ph-^lp'h'- 1 ). (1) 

After this replacement, the numerical solution of the Continuum RPA equations follow the 
lines as outlined in Ref. ||14j| . The orthogonalization of the particle states on the hole states 
is performed by Schmidt's procedure. 

With the exception of the T = 1 triad of bound states in the A = 12 nuclei, all other 
final states (which contribute to the inclusive cross section) lie above the particle emission 
threshold and therefore are described as resonances within the continuum RPA model. As 



residual interactions we used as before the finite-range G-matrix |T5j derived from the Bonn 
NN potential (BP) and, in order to estimate the theoretical uncertainties, also an empirical 
zero-range Landau-Migdal force (LM). 

To derive an estimate of the partial occupancies of the pi/ 2 subshell, we have performed 
a shell model calculation for 12 C with the OXBASH code [|16l and the Cohen-Kurath inter- 



action ||17|| . This calculation suggests a partial occupation for the pi/ 2 subshell of 0.75 for 



protons and neutron, respectively. This result is rather close to the value 0.8 used in Ref. 



11 



First we recalculated the various exclusive processes leading to the T = 1 triad in the 
A = 12 nuclei. Since these states are discrete we used the standard RPA in this calculation, 
again including partial occupancies as outlined above. As single particle states we considered 
the complete (ls)(lp)(2sld)(2plf) model space. For the proton and neutron p 3 / 2 and p\/ 2 
single particle energies we used the experimental values; the others have been determined 
from a Woods-Saxon potential which describes the 12 C charge distribution well 18|| . We 



have checked that an enlargement of the model space does not affect our results. Our RPA 
results are compared to data in Tables 1 and 2. 

In Table 1 the rescaling factors N needed to bring calculation and data in agreement 
are shown. Note, that allowing the partial occupancy of the pi/ 2 subshell reduces the value 
of N substantially, from about 4 to less than 2. At the same time the rescaling factors are 
essentially identical for the three processes considered, showing consistency. The necessity 
of using such renormalization reflects again the well known fact that random phase approxi- 
mation calculations do not describe all (proton-neutron) correlations necessary to reproduce 
the quenching of the Gamow- Teller strength This aim can be only achieved in 
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shell-model calculations within a complete major shell, introducing additionally the renor- 
malization of the spin operator by the universal factor 0.77 [28- |30|l . Deviating from this 



universal finding, the Cohen-Kurath interaction had been designed to describe the weak 
interaction transition rates in 12 C without incorporating the renormalization of the spin 
operator. Thus the fact that iV > 1 in our calculation again shows the wellknown result 
that the reproduction of the Gamow- Teller strength requires the consideration of genuine 
correlations beyond a simple mean-field approach. 

In Table 2 we compare the results on the exclusive neutrino induced reactions with 
data. The entries are based on J\f new = 1.5 and N old = 4.0. The agreement remains good, 
and at the same time, as mentioned above, the required renormalization is substantially 
smaller than in our previous RPA calculation. The reduction is caused by the destructive 
interference of the P3/2 and pi/2 configurations in the RPA description of the Gamow- Teller 



transition, as discussed by Auerbach et al. |TT 



We then studied the three inclusive processes within our improved continuum RPA 
model. The noticeable mixture of pi/2 configuration into the 12 C ground state allows for the 
excitation of a larger number of resonances than in the previous calculation. This is demon- 
strated in Fig. D, where the 12 C(^ At , ^t~) 12 N* cross section (calculated with the BP force) is 
shown as a function of the excitation energy uj of the nucleus. However, this figure also 
demonstrates that the strength going into these additional resonances is taken away from 
those transitions already present for the pure p^/2 ground state. In fact, Table 3 shows that 
the total inclusive rate and cross sections are only slightly reduced by the mixing of a pi/2 
component into the ground state. The inclusive muon capture rate and the 12 C(z/ e , e~) 12 N* 
cross section agree very well with data: For muon capture the improved continuum RPA 
yields for both residual interactions rates close to the experimental value (32.8 • 10 3 s _1 ), 
and the calculated 12 C(u e , e~) 12 N* cross sections are also close to the most recent data 
(5.1 ± 0.6 ± 0.5) ■ 10~ 42 cm 2 § and (5.7 ± 0.6 ± 0.6) ■ 10" 42 cm 2 @. Our calculation of 
the inclusive 12 C(^, /i~) 12 N* cross section results also in a reduced value compared to our 
previous study, where the reduction is about half (0.8 ■ 10~ 40 cm 2 for BP, 1.9 ■ 10~ 40 cm 2 
for LM) due to the revised flux and half (0.7 ■ 10" 40 cm 2 for BP, 0.9 • 10~ 40 cm 2 for LM) 
due to the partial occupation of the pi/ 2 subshell in the ground state wave function. Our 
new total 12 C(i^,, /z~) 12 N cross section (inclusive + normalized exclusive) is then 17.8 ■ 10~ 40 
cm 2 and 17.5 • 10~ 40 cm 2 for the BP- and LM-force, respectively, to be compared with the 
last experimental value of ((12.4 ± 0.3(stat.) ± 1.8(syst.) • 10~ 40 cm 2 ). Thus the disturbing 
discrepancy is less severe than before, but it not entirely removed. 

Our result does not support in detail the argumentation in Ref. |TT| that the discrepancy 
between RPA calculations and data can be essentially fully removed if pi/2 partial occupancy 
is considered (Ref. |Tl| gives (13.5 — 15.2)-10~ 40 cm 2 ) and all inclusive and exclusive processes 



can be simultaneously and consistently described. We also mention that Singh et al. |32 
have studied the various processes based on the local density approximation and find a total 
(exclusive + inclusive) muon capture rate of (3.6 ± 0.22) ■ 10 4 s~ x (the experimental 



value for the total rate is (3.8 ± 0.10) • 10 s ) and a C(z/ M , fM~) N cross section of (16.5 ± 
1.3)10 -40 cm 2 ||32|| , slightly lower, but consistent with our new result. Notably, Mintz and 
collaborators predicted the 12 C(f M , /i~) 12 N cross section |54| , but the model assumptions 
within the elementary particle model used in the approach of these authors are questionable 
and not justified for this reaction ]T2|. 
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In summary, we have found that the consideration of a partial pi/2 occupancy lowers the 
exclusive cross sections and muon capture rate noticeably, but it reduces the results for the 
inclusive processes only slightly. The reason for the reduction of the exclusive cross sections 
has been discussed above. But can one understand why the continuum RPA is apparently 
able to describe the inclusive processes without the need for parameter adjustment, unlike 
the exclusive Gamow- Teller transitions? Another way of thinking about this problem is the 
question to what extent the correlations of nucleons within the p-shell influence the results. 
To shed the light on this, one can evaluate the total strength of various operators O, i.e. the 
norm of the vector 0\g.s.), with and without the effect of correlations. Note that the total 
strength depends only on the ground state wave function and is therefore relatively easy to 
evaluate. (Here we follow the points made in Ref. [|3~5|D . 

For the positive parity operators this is done in Table 4. Table 4 illustrates the well 
known fact that the p shell correlations are very important for the Gamow- Teller operator 
ar. (Note that the RPA gives the total Gamow- Teller strength of 5.5, only slightly reduced 
when compared to the 'naive' estimate of the independent particle model. Also, the exact 
shell model predicts that a strength of 0.12 goes to excited 1 + states, which are absent in 
the naive model.) 

But the situation with the quadrupole operators is rather different. The total p shell 
strength of the spin-independent operator, and the strength summed over the multipoles 
of the spin-dependent operators is affected by the correlations only on the 30 — 40% level, 
even though the individual spin-dependent multipoles are affected more. Moreover, the p 
shell strength represents only a small fraction of the total quadrupole strength, which is 
concentrated in the 2hu excitations, unaffected by correlations as long as we assume that 
the ground state has only p shell nucleons. 

The lowest 2 + resonance has a dominating (Op) configuration and thus its contribution 
to the various cross sections might be affected by the present more flexible description of the 
12 C ground state. We find that the contribution of the 2 + multipole to the z/ e -induced cross 
section and to the muon capture rate increases slightly (by 0.02 ■ 10~ 42 cm 2 and 0.17 • 10 3 
s _1 , respectively), but is still negligible compared to the total cross section and rate (0.8% 
and 2.8%, respectively). For the z^-induced reaction, the effects of the partial occupancy 
formalism and of the improved experimental neutrino flux nearly cancel each other so that 
the partial cross section of the 2 + multipole slightly decreases (by 0.08 • 10 -40 cm 2 ). As the 
total cross section is also reduced, the contribution of the 2 + multipole, of which the (Op) 
configuration is only a small fraction, to the total cross section is increased to 12% in our 
present calculation. 

However, the inclusive reactions we are considering are dominated by the excitations 
of the negative parity states. To what extent does the strength depend on the occupation 
of the pi/2 subshell for these operators? According to Ref. [I]J the cross sections should 
be noticeably reduced when there are about 1.6 nucleons in the pi/2 subshell. In order 
to test this point, we plot in Fig. || the dipole and quadrupole strength as a function of 
the pi/2 subshell occupation. We find that when summed over multipoles the strength is 
totally independent of this occupation number. But even the individual multipoles depend 
on the occupation numbers only mildly. We find that the p shell correlations are relatively 
unimportant for the full strength. 

Of course, the partial occupation of the pi/ 2 subshell leads to the occurrence of new 
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states, reshifting the strength. As the pi/ 2 orbital reflects an excited state compared to 
the j>3/2 orbital, the strength is shifted slightly to higher energies. Thus, even if the total 
strength is unchanged, this redistribution of strength changes the cross section due to the 
related change of kinematics. As the strength is shifted to higher energies, its weight in 
the cross section is reduced, as apparent from our results. The relative importance of this 
effect decreases with increasing momentum and energy transfer, explaining why we find a 
reduction of about 20% for the v e induced experiment, while it is only 4% for the ^-induced 
reaction. 

The LSND collaboration has recently conducted a search for — > v e oscillations for 
high energetic (E u < 300 MeV) stemming from 7r + decay in flight. The oscillation signal 
consists of isolated, high-energy electrons (60 MeV< E e - <200 MeV) in the detector coming 
from charged current 12 C(v e , e~) 12 N reactions. Such events cannot be caused by zvneutrinos 
from pion decay at rest, which are also produced in the beam stop at LAMPF, because their 
energies is lower than 52.8 MeV. In the experiment an excess in the number of these events 
was found, which, if interpreted as an oscillation signal, leads to an oscillation probability 
of (2.6 ± 1.0 ± 0.5) x 10 -3 p6|| . This value and its systematic error depends dominantly on 
the knowledge of the inclusive 12 C(is e , e~) 12 N cross section and the absolute neutrino flux 
through the detector. To aid the analysis, we present in Table 5 this cross section calculated 
within the improved RPA model. Comparing with our previous results in Ref. [[Hj, we find 
for most of the neutrino energies a reduction of the cross section between 5-10%. As this 
reduction is smaller than the 10%-error attributed to the uncertainty of 12 C(u e , e~) 12 N cross 
section in the LSND-analysis it does not require a re-analysis of the experiment. 

Acknowledgments We were supported in part by the Swiss National Science Foundation, 
the Danish Research Council and by the U.S. Department of Energy under grant DE-FG03- 
88ER-40397. 
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TABLES 



TABLE I, The rescaling factors N of the matrix elements for the exclusive transitions were 
determined for both of the applied residual interactions by comparison of the calculated rates 
(c^pp", to^yi) with the experimental beta-decay rates and partial muon capture rate to the 12 B 
ground state (in units of s" 1 ). Also shown are the previous results (u; B p, w£m) obtained for a 
completely filled p 3 / 2 subshell. 



process 


data Ref. 


. .new , .new t\t2 at2 
^BP ^LM iV BP LM 


, .old , .old aj-2 nj2 
W BP W LM iV BP LM 


12 C(/i-,^) 12 B g . s . 

12 B g . s .(r) 12 C g . s . 
12 N g .,(/?+) 12 C g . s . 


6 050±300(*) f5|- 
33.36 ±0.13 
59.58 ±0.46 


22] 
23] 
23] 


9 330 8 538 1.54 1.41 

49.6 44.3 1.49 1.33 

98.7 88.2 1.66 1.48 


22 780 25 400 3.77 4.20 
123.6 128.8 3.71 3.86 
247.1 257.4 4.15 4.32 



TABLE II. The exclusive cross sections for charged and neutral current neutrino scattering (in 
units of 10 -42 cm 2 ). The results of our improved RPA calculation for both of the applied residual 
interactions (dpp*, o'lmO are compared to the data and the previous RPA calculation (<7 B p, o"lm) 
without partial pi/ 2 subshell occupation. 



process 


data Ref. 


_new _new 
a BP "LM 


old old 
°BP "LM 


12 C(z, e ,e-) 12 N g . s . 
12 C( W ") 12 N g . s . 


10.5 ± 1.0(stat.) ± 1.0(syst.) || 
8.9 ± 0.6(stat.) ± 0.75(syst.) § 
9.1±0.4(stat.)±0.9(syst.) § 
66 ± 10(stat.) ± 10(syst.) § 


8.9 8.9 

68 73 


9.3 9.3 

63 63 


12 C(i/, i/) 12 C*(15.11) 
12 C(^) 12 C*(15.11) 


10.4 ± 1.0(stat.) ±0.9(syst.) 
3.2±0.5(stat.) ±0.4(syst.) \'. 


§ 

15] 


10.5 10.5 
2.8 2.7 


10.5 10.6 
2.8 2.8 



TABLE III. The inclusive muon capture rate to, (in 10 3 s _1 ) and the cross sections a for the 
12 C(v e , e~) 12 N* (in units of 10 -42 cm 2 ) and the total (inclusive + exclusive) cross section for the 
12 C(z^j, /U~) 12 N (in 10 -40 cm 2 ) reactions.. The results of our improved continuum RPA calculation 
((uj/a) new ) are compared to the data and the previous ((w/cT) old ) continuum RPA calculation 



without partial p 1 / 2 subshell occupation. 



process 


data Ref. 






12 C(^",^) 12 B* 
12 CK,e-) 12 N* 

12 C( W -) 12 N 


32.8 ±0.8 H 
5.1 ±0.6 ±0.5 § 
5.7 ±0.6 ±0.6 § 
12.4 ±0.3 ±1.8 [13] 


32.7 31.3 
5.4 5.6 

17.8 17.5 


34.2 33.3 
6.3 5.9 

19.3 20.3 



7 



TABLE IV. The full strength within the nuclear p shell evaluated for the operators in column 
1. The 'SM' column is the exact shell model result calculated with the Cohen-Kurath interaction. 
The column named 'naive' corresponds to the (P3/ 2 ) S configuration. In the last column the strength 
for transitions with 2hu) is shown for comparison. 



Operator 


SM 


naive 


2hu 


Gamow- Teller (err) 


1.51 


8.00 


0.00 


r 2 Y 2 


1.37 


1.98 


9.95 


r 2 (Y 2 a) I=l 


0.11 


0.08 




r 2 (Y 2 aY= 2 


0.33 


0.75 




r 2 (Y 2 aY= 3 


0.20 


0.00 




ExT 2 (Y 2 a^ 


0.64 


0.83 


7.47 
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TABLE V. Total 12 C(z/ e , e") 12 N, 12 C(P e , e+) 12 B, and exclusive cross sections to the 12 N and 
12 B ground states for a mesh of neutrino energies E Ve . The cross sections have been calculated 
with the BP-interaction and are given in units of 10~ 42 cm 2 . Energies are in MeV, exponents are 
given in parentheses. 



K e 


12 C(z/ e ,e-) 12 N 


12 C(z/ e ,e-) 12 N 9S 


12 C(P e ,e+) 12 B 


12 C(v e ,e+) 12 B gs 


20 


2.85 (-1) 


2.84 (-1) 


8.00 (-1) 


7.90 (-1) 


30 


5.65 


4.90 


6.05 


5.01 


40 


2.23 (+1) 


1.46 (+1) 


1.84 (+1) 


1.17 (+1) 


50 


5.99 (+1) 


2.79 (+1) 


4.14 (+1) 


1.95 (+1) 


60 


1.31 (+2) 


4.31 (+1) 


7.80 (+1) 


2.72 (+1) 


70 


2.47 (+2) 


5.82 (+1) 


1.30 (+2) 


3.41 (+1) 


80 


4.21 (+2) 


7.16 (+1) 


1.99 (+2) 


3.97 (+1) 


90 


6.63 (+2) 


8.25 (+1) 


2.84 (+2) 


4.39 (+1) 


100 


9.79 (+2) 


9.03 (+1) 


3.86 (+2) 


4.69 (+1) 


110 


1.37 (+3) 


9.49 (+1) 


5.01 (+2) 


4.90 (+1) 


120 


1.85 (+3) 


9.70 (+1) 


6.29 (+2) 


5.04 (+1) 


130 


2.41 (+3) 


9.70 (+1) 


7.68 (+2) 


5.14 (+1) 


140 


3.06 (+3) 


9.58 (+1) 


9.17 (+2) 


5.22 (+1) 


150 


3.79 (+3) 


9.37 (+1) 


1.08 (+3) 


5.29 (+1) 


160 


4.61 (+3) 


9.15 (+1) 


1.24 (+3) 


5.35 (+1) 


170 


5.51 (+3) 


8.94 (+1) 


1.42 (+3) 


5.41 (+1) 


180 


6.49 (+3) 


8.77 (+1) 


1.59 (+3) 


5.47 (+1) 


190 


7.54 (+3) 


8.63 (+1) 


1.78 (+3) 


5.52 (+1) 


200 


8.66 (+3) 


8.51 (+1) 


1.96 (+3) 


5.58 (+1) 


210 


9.84 (+3) 


8.42 (+1) 


2.15 (+3) 


5.63 (+1) 


220 


1.11 (+4) 


8.35 (+1) 


2.34 (+3) 


5.68 (+1) 


230 


1.23 (+4) 


8.28 (+1) 


2.53 (+3) 


5.72 (+1) 


240 


1.36 (+4) 


8.22 (+1) 


2.72 (+3) 


5.77 (+1) 


250 


1.50 (+4) 


8.17 (+1) 


2.91 (+3) 


5.81 (+1) 


260 


1.63 (+4) 


8.11 (+1) 


3.11 (+3) 


5.85 (+1) 


270 


1.76 (+4) 


8.06 (+1) 


3.30 (+3) 


5.89 (+1) 


280 


1.90 (+4) 


8.01 (+1) 


3.50 (+3) 


5.93 (+1) 


290 


2.03 (+4) 


7.96 (+1) 


3.70 (+3) 


5.97 (+1) 


300 


2.16 (+4) 


7.91 (+1) 


3.90 (+3) 


6.00 (+1) 
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FIGURES 



FIG. 1. The 12 C(f^, / u~) 12 N* cross sections as a function of the excitation energy u> of the 
nucleus. The present results obtained with partial occupancy of the pi/ 2 orbital (labelled 'new') 
are compared to those assuming a pure (P3/2) 8 closed-shell configuration for the 12 C ground state 
(labelled 'old'). 

FIG. 2. Total strength of the dipole (upper part) and quadrupole (lower part) operators versus 
the occupation of the pi/ 2 subshell. The curves are labeled by the corresponding angular momentum 
of the r l (Yi<j) I=x operators for / = 1 and I = 2. 
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